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SUMMARY 


A simplified fourwall interference assessment method has 
been described, and a computer program developed to facilitate 
correction of the airfoil data obtained in the Langley 0.3-m 
Transonic Cryogenic Tunnel (TCT) . The procedure adopted is 
to first apply a blockage correction due to sidewall boundary- 
layer effects by various methods. The sidewall boundary- 
layer corrected data are then used to calculate the top and 
bottom wall interference effects by the method of Capallier, 
Chevallier and Bouinol, using the measured wall pressure 
distribution and the model force coefficients. The inter- 
ference corrections obtained by the present method have been 
compared with other methods and found to give good agreement 
for the experimental data obtained in the TCT with slotted top 
and bottom walls. 


INTRODUCTION 

Wall interference is a problem of concern in testing of 
airfoils in wind tunnels, particularly at transonic speeds. 
This is largely due to the complex flow features at the 
ventilated walls introducing uncertainties in the boundary 
conditions. Therefore, modern wall interference calculation 
methods, in lieu of the classical boundary conditions, use 
experimentally measured values of the pressure and/or flow 
inclinations at the wall to correct the test data. 

A comprehensive review of the different methods of 
calculating the two-dimensional wall interference has been 
given recently by Mokry et al (ref. 1) . Table I lists the 
salient features of some of these methods. Most of them are 
based on subsonic flow theory and still give useful results in 
the low transonic regime as long as the flow is subcritical at 
the walls. The only methods which employ transonic analysis 
are those of Murman (ref. 2) and Kemp (ref. 3) . 
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Kemp's method employs numerical solution of the transonic 
small disturbance equation, and requires as input the measured 
pressure distributions on the airfoil model, and near the top 
and bottom walls. The solution of this inverse problem gives 
an effective shape of the airfoil, which accounts indirectly 
for the viscous effects also, on the airfoil. This effective 
airfoil shape is then used to calculate the unbounded or free 
air solution. The Mach number and the angle of attack are 
iterated to give the best match between the calculated and the 
measured pressure distributions. 

Kemp's method can be considered to represent the state of 
the art in interference calculations and can be used as a 
standard of comparison against which other simpler methods can 
be validated. 

Another method of analysis which uses the panel technique 
to calculate numerically the entire interference velocity 
field is due to Smith (ref. 4) . 

Most of the literature on interference calculation in 
two-dimensional wind tunnels deal with top and bottom wall 
influence which is inviscid in nature and is uniform across 
the span. It is generally assumed that the flow in the wind 
tunnel is nearly two-dimensional. Another source of inter- 
ference which has been receiving considerable attention 
recently is the influence of the sidewall boundary-layers 
which is essentially viscous in nature and introduces three- 
dimensional perturbations across the width of the tunnel. 

Though, it was recognised by Preston (ref. 5) in 1944 
that the interaction of the sidewall boundary-layers with the 
airfoil flowfield can cause departure from two-dimensional 
conditions, only recently there has been systematic effort 
towards understanding the extent of three dimensional 
influence and means to correct for the same. Most noticeable 
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development has been due to Barnwell (ref. 6) who considered 
the changes in the sidewall boundary-layer to introduce 
crossflow velocities across the width of the test-section. 
Assuming linear variation of the crossflow velocity, and using 
a simplified treatment of the sidewall boundary-layer growth, 
Barnwell suggested a simple correction in the form of a 
modified Prandtl-Glauret rule, in terms of the test section 
Mach number and the undisturbed values of the sidewall 
boundary-layer displacement thickness and shape parameter. 
Sewall (ref. 7) extended this approach to include transonic 
effects by using the von Karman similarity rule. A modified 
form of the correction was proposed recently by the present 
author (ref. 8) , by considering the sidewall boundary-layer 
to cause changes in both the Mach number and the airfoil 
thickness. Also, a new approach which accounts for the 
nonlinear variation of the crossflow velocity across the 
width of the tunnel has been suggested and the correction to 
the test Mach number has been shown to be a function of the 
airfoil aspect ratio (ref. 9) . 

A four-wall interference assessment procedure in two- 
dimensional airfoil testing was suggested by Kemp and Adcock 
(ref. 10) by combining the Barnwell -Sewall sidewall boundary- 
layer correction with the Kemp's method for top and bottom 
wall interference calculations. This four-wall interference 
assessment code (ref. 11) has been used extensively to make 
post-test corrections for the airfoil test data obtained in 
the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel (TCT) . 
Since, this method requires significant computational effort, 
its use is recommended for selected data points requiring 
accurate treatment. 

For routine calculation of interference to all the data 
points in an airfoil test, it is desirable that a simpler 
approach amenable for quick calculation is employed. 
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Particularly, when the interest is confined to making global 
corrections in the region of the airfoil, the methods due to 
Capalier et al. (ref. 12), Mokry (ref. 13), Sawada (refs. 14, 
15) , and of Ashill and Weeks (refs. 16, 17) are attractive. 
These methods require little computational effort and are well 
suited for on-line calculations, if necessary. 

The method of Ashill and Weeks requires measured pressure 
distribution and flow inclination at the walls, and is well 
suited for solid wall tunnels where flow inclinations at the 
wall are known more accurately. The corresponding Cauchy 
Integral formula is solved, and no model representation is 
required. 

For ventilated wall tunnels, where the flow inclinations 
are difficult to measure, the other methods are convenient. 
These methods use the measured pressure distribution near the 
wall in conjunction with the model represention by appropriate 
singularities. Mokry' s method solves the potential flow 
problem in the rectangular domain bounded by the testsection 
geometry. The methods of Capalier et al . , and of Sawada, 
though based on different approaches, solve the problem in the 
infinite strip bounded between the top and bottom walls, and 
give similar results to a first approximation. All these 
methods, including the Kemp's method have been found to give 
nearly identical interference corrections for the BGK-1 
airfoil test data obtained in the NAE 15"x60" tunnel with 
perforated top and bottom walls (ref. 1) . 

The good agreement between the various methods for the 
NAE test case prompted the present investigation to examine 
the application of the method of Capalier et al., to the 
airfoil test data obtained in the 0.3-m TCT two-dimensional 
test section with slotted top and bottom walls. In addition 
to facilitate calculation of fourwall interference effects, 
the sidewall boundary-layer correction methods of references 
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6-9, were considered in a sequential manner. The approach 
was validated by comparing results obtained with Kemp's method 
for some test cases, and good agreement was observed. This 
suggests that for routine calculation of interference correc- 
tions in the 0.3-m TCT airfoil tests, the method of Capalier 
et al., can be successfully employed. This method has been 
implemented in a computer program to facilitate calculation of 
interference corrections directly from the airfoil test data 
files, either for sidewall or top/bottom wall, or both. The 
objective of this report is to present the method of analysis 
used, its validation results and details of the computer 
program developed for correcting the 0.3-m TCT test data. 

NOMENCLATURE 


A 

b 



M c 

M, 

M, 

M+ 


c,b 

c,m 


^c 

Re 


m 


*w' 


'w 


non-dimensional area of the model 

width of the tunnel 

model chord length 

drag coefficient 

lift coefficient 

pressure coefficient 

functions defined by eqns. (12) and (13) 

total height of the tunnel 

sidewall boundary-layer shape factor 

= (2§*/b) (2 +1/H -M t 2 ) 

modified form of k (see equation 6) 

= (tt/ 3 b/i) 

wave length in terms of model chord 
corrected Mach number 

corrected Mach number by Barnwell-Sewall Method 
corrected Mach number by Murthy's method 
test Mach number 

Reynolds number based on model chord 
freestream Reynolds number/meter 
non-dimensionalised interference velocities in 
x and y direction due to top and bottom walls 
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u 


m 


v i 

Wi, w 2 

Y 


x, y 
a 

P 

7 

8 * 

i 

M 

a 

AM 

Aa 


non-dimensional ised induced velcoity due to model 

tunnel upstream flow inclination 

Functions defined by equations (18) and (21) 

= (3h/2 

streamwise and normal coordinates 
angle of attack 
= (1 - M 2 ) 1 / 2 
vortex strength 

displacement thickness (sidewall boundary-layer) 

velocity potential 

doublet strength 

source strength 

correction to Mach number 

correction to angle of attack 


Subscripts: 


e refers to end of the testsection 

s refers to beginning of the testsection 

1 testsection lower wall 

u testsection upper wall 

w refers to wake 

DESCRIPTION OF THE METHOD 
Sidewall Boundary-Layer Correction 

The method adopted in the present analysis to calculate 
the fourwall interference effects is to apply first the side- 
wall boundary-layer correction. The option to use either the 
method of Barnwell-Sewall , or Murthy has been considered. 
Briefly both the methods employ the simplified assumption that 
the sidewall boundary-layer growth can be represented by von 
Harman' s momentum integral equation and the effect of skin 
friction can be ignored in comparison with the model induced 
pressure gradients. 
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In the Barnwell-Sewall method, it is further assumed that 
the sidewall boundary-layers induce spanwise velocities which 
vary linearly across the width of the test-section. Using the 
von Kantian' s similarity parameter, a correction to the 
measured test-section Mach number is suggested in terms of the 
empty tunnel sidewall boundary-layer displacement thickness 
and the shape factor. The value of the corrrected Mach number 
M c b is given by (ref. 7) 

(1 - M 2 cb ) 3/4 (1 - M 2 t + k ) 3 / 4 

M^b = »t U) 

where M t corresponds to the test Mach number, and k is a 
constant calculated using the empty tunnel boundary-layer 
displacement thickness ( 5 *) and the shape parameter H. 

k = (2 5*/b) (2 + 1/H -M 2 t ) ( 2 ) 

The measured pressure and force coefficients are multiplied by 
the factor (ref. 8 ) 


= (M t /M C/b ) 2/3 (3) 

to give the corrected values. 

It may be noted that in the Barnwell-Sewall method, the 
correction to the test Mach number is proposed only for tran- 
sonic speeds. In a form proposed earlier by Barnwell (ref. 

18) , the correction to the test Mach number was not defined at 
lower Mach numbers. This deficiency was overcome by Murthy 
(ref. 8 ) , by using a modified coordinate transformation of the 
govering small disturbance equation. It was shown that the 
flow in the wind tunnel with sidewall boundary-layers can be 
considered as an equivalent two-dimensional flow over an 
airfoil of reduced thickness at a reduced Mach number. With 


7 


this approach, the corrected Mach number M c m can be expressed 
in a simplified form as (ref. 8) 

M c,m = “t/* 1 + k ) V2 ( 4 ) 

The corrected Mach number M c , m given by equation (4) is 
valid from subsonic speeds to transonic speeds, and agrees 
with the Barnwell-Sewall correction at transonic speeds for 
small values of k. The corresponding correction factors for 
the measured pressure and force coefficients are 


for 

subsonic speeds : 

= (1 

+ k) 1 / 2 

(5a) 

for 

transonic speeds: 

= (1 

+ k) V3 

(5b) 


The difference between the corrections obtained by the 
Barnwell-Sewall method and the Murthy's method is not 
significant at transonic speeds, and either of them can be 
considered to be equally valid within the small disturbance 
approximations made in deriving the corrections. The method 
of reference (8) , however, appears much simpler and facilita- 
tes a continuous correction from subsonic to transonic speeds. 

Aspect Ratio Effects 

Methods discussed above for calculating the sidewall 
boundary-layer effects are not entirely satisfactory. The 
corrections derived depend only on the sidewall boundary-layer 
parameters, model span and the test Mach number, and are 
independent of the model chord. They can be considered to be 
applicable only when the reduced aspect ratio (= (3 b/c) of the 
model is small, so that the effect of the sidewall boundary- 
layers is nearly one-dimensional, at least in the vicinity of 
the airfoil. For higher aspect ratio models, one can expect 
the effects to be much smaller in the mid-span region of the 
model where the measurements are made. To represent this 
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diminishing effect of the sidewall boundary-layers with 
increasing aspect ratio, an improved correction suggested by 
the present author (ref. 9), has been used. This method, 
instead of the linear crossflow velocity assumption, considers 
the two-dimensional flow between a wavy wall and a fixed wall 
to represent the cross flow velocity effects. With this wavy 
wall approach, the effect of aspect ratio is included by 
defining the constant k in equations (1) to (5) , in a modified 
form given by 

k m = (2 5*/b) (2 + 1/H - M 2 t ) (k 2 /Sinh(k 2 ) } (6) 

where k 2 = ( n (5 h/l ) (7) 

The factor k 2 /Sinh(k 2 ) depends on the test Mach number, 
the model span and a length scale i representing the model 
chord. In the limit of vanishing aspect ratio, this factor 
approaches a value of one, and the equation (6) reduces to 
equation (2) . 

When considering the aspect ratio correction, it is 
necessary to define what constitutes a typical length scale 
in terms of the model chord. This has been examined in 
reference (9) , and it appears that a value of i =2c appears 
reasonable considering the fact that the effect of the airfoil 
on the sidewall boundary-layers is distributed over a distance 
of about twice the chord of the airfoil . It may be noted that 
this aspect ratio correction is based on two-dimensional 
analysis. Hence, it is likely that the reduction in the 
sidewall boundary-layer effects can be smaller than predicted 
due to three-dimensional nature of the flow at the junction. 

It may be noted that the sidewall boundary-layer correc- 
tion methods discussed account only for the negative blockage 
effect caused by the thinning of the boundary-layer in the 
airfoil region due to favourable pressure gradient. It is 
likely that downwash effects can be present as evidenced by 
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detailed measurements on a Cast-7 airfoil over a wide range of 
aspect ratios (ref. 19) . However, the physical mechanism 
causing such effects does not appear to be well understood and 
cannot be represented by simple mathematical models. A 
detailed discussion of the physical phenomenon associated with 
the sidewall boundary-layer effects has been presented by 
Winter and Smith (ref. 20), and it now appears that the 
changes in the boundary-layer displacement thickness suggested 
by Barnwell appears the most plausible one. Considering the 
uncertainties in angle of attack in two-dimensional airfoil 
testing, the simple sidewall boundary-layer correction is 
quite useful, particularly in ventilated wall tunnels designed 
for minimum blockage (ref. 21) . 

Sidewall Boundary-Layer Parameters 

To apply the sidewall boundary-layer correction, it is 
necessary to know the empty tunnel boundary-layer displacement 
thickness and the shape factor at the location of the model. 
These parameters are generally measured during tunnel calib- 
ration. If the measured values are not available, theoretical 
estimates can be made by assuming a fictitious flat plate 
boundary-layer growth (ref. 22). For the 0.3-m TCT, measure- 
ment of the sidewall boundary-layers have been made at various 
times using a wall mounted rake located upstream of the model 
station (see ref. 23) . Using these measurements, the values 
at the model station were estimated using flat boundary-layer 
theory. Based on these calculations, the following empirical 
formulae have been used in the program for calculation of 
boundary-layer parameters. 

5* (mm) = 6.42266 - 0.59613 log(Re m ) + 

M t {0. 44608 -0.01333log(Re m ) } (8) 

H =1.54608 +0.44299M t -0 . 04828log (Re m ) (9) 
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TOP AND BOTTOM WALL INTERFERENCE 


The method used in the present report to calculate the 
top and bottom wall effects is that due to Capalier et al. 

The details of the method are given in references (1) and 
(12) . In this method, the pressure coefficients at or near 
the top and bottom walls from upstream infinity to downstream 
infinity, are used to prescribe the boundary conditions to 
solve the interference problem. For an airfoil located 
midway between the two walls of a tunnel (fig. 1) , the 
components of the non-dimensionalised interference velocities 
u w and v w , in the x and y directions at the location of the 
model are given by 


OO 

u w (0,0) = (1/ |3h)J’[f u (x) + f x (x) ]/[2cosh( 7r x/0h) ] dx (10) 

— OO 
OO 

v w (0,0) = v i +(l/h)/(f u (x) - f 1 (x)]/(l+exp(27rx//3h)] dx (11) 

— OO 

where f u and f^ are functions of the measured pressure coeffi- 
cients Cp, u and Cp^ on the upper and lower walls, and the 
non-dimensionalised induced velocities u m u and u m ^ due to 
airfoil model in free air. The corresponding expressions for 
f u and fj are 


f u = -(l/2)C p/U (x) 

" u m,u< x ) 

(12) 

f x = -(l/2)C p ^ 1 (x) 

- V:l(x) 

(13) 


The term v^ in equation (11) is a constant, and refers to 
the tunnel upstream flow inclination which has to be arrived 
at by empty tunnel calibration. The induced velocities u m u 
and u m near upper and lower walls due to airfoil model in 
free air are obtained by representing the model by appropriate 
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singularities: doublet for model blockage, source for wake 
blockage and vortex for lift effects. The corresponding 
velocity potentials are given by 


doublet: 

0 

= (]U/27r(3) [x/(x 2 +0 2 y 2 ) ] 

(14a) 

source : 

0 

= (0/4tt(3) log (x 2 +0 2 y 2 ) 

(14b) 

vortex : 

0 

= -(7/2 K ) tan -1 (py/x) 

(14c) 


By differentiating equations 14a-c, the corresponding 
model induced velocities at the top and bottom wall locations 
(y=±h/2) can be written as 


doublet 

: u m,u= 

u m , 1= (' -M/ir 0) [ ( x 2 - Y 2 ) / ( x 2 +Y 2 ) 2 ] 

(15a) 

source 

: u m,u= 

u ro , l = (°/ 2 7Tj3) [X/(X 2 +Y 2 )] 

(15b) 

vortex 

: u m , u = 

” u m , i = (7/2 7T ) [Y/(x 2 +Y 2 )] 

(15c) 

where 

Y= /3 h/2 . 



For a thin airfoil at incidence, the strength of the 
various singularities are determined by the model area of 
cross-section, and the lift and drag force coefficients. 


doublet 

: M 

= A c 2 

(16a) 

source 

: a 

(N 

0 

0 

II 

(16b) 

vortex 

: y 

= (C c ± )/2 

(16c) 


The calculation of the Mach number and the angle of 
attack corrections is quite straightforward provided the 
measured wall pressure data are available over a sufficiently 
long distance both upstream and downstream of the model, so 
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that the errors in integration due to truncation is small. 
However, this is a factor that is largely governed by the 
available length of the tunnel test-section, and suitable 
extrapolations may have to be made. 

Correction for Mach number 

Substituting equations (12) and (13) , the equation (10) 
for the blockage interference velocity at the model u w (0,0) 
can be written as 

u w (0,0)= -(1/ h) J[(C pfU +Cp fl )/2 +(u m/U +u mfl )]W 1 dx (17) 

- OO 

where Wj(x) = l/[ 2 cosh( 7 rx//?h) ] (18) 

The function W-^(x) can be considered as a weighting 
factor which multiplies the measured wall pressures, and the 
model induced velocities. The variation of W-^fx) along the 
length of the testsection is shown in figure (2) , for Mach 
numbers 0 . 0 and 0.7. Its value is maximum at the center of 
the testsection and decays exponentially both upstream and 
downstream for large (x/h) . Beyond x/h = ±2, its contribu- 
tion is negligible, and the blockage correction will be 
insensitive to uncertainties or errors in wall pressure 
measurements in this region. Also, with increase in the test- 
section test Mach number, the contribution to correction is 
limited to a narrower region about the test-section center- 
line. Hence, the exponential behaviour of the weighting func- 
tion W-l has the beneficial effect that the integration of 
equation (17) can be limited between the upstream and 
downstream ends of the test-section without loss of accuracy. 
The test Mach number corrected for the top and wall inter- 
ference is given by (ref. 1) 

M c = M t [l + (1 + 0.2 M t 2 )u w (0,0)] (19a) 
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The force coefficients are referenced to the corrected Mach 
number M c , by multiplying the measured values by the factor 


M t 2 [1 + 0.2 M c 2 ] 


M c 2 [1 + 0.2 M t 2 ] 


(19b) 


Correction for Angle of Attack 

Substituting equations (12) and (13), the equation (11) 
for the interference velocity v w (0,0), at the model can be 
written as 


v w (0,0)= -d/h)/r<C p , u -C pil )/2 +<u miU -u m;1 )]W 2 dx (20) 

— oo 

where 


W 2 (x) = 1/[1 +exp(2 7Tx/j3h) ] (21) 

The variation of the weighting function W 2 (x) with x is 
shown in figure (3) , and it may be noted that it has an 
asymptotic value of one upstream and zero downstream. The 
variation between these two limits occurs over a narrow 
region extending about one testsection height either side. 

This suggests that the integration of equation (20) can be 
truncated at a suitable location downstream, since the 
contribution beyond that region will be negligible. 

However, the same argument does not apply for the upstream 
end, since the weighting function W 2 (x) is almost equal to one 
beyond about one testsection height. Hence, the contribution 
to the integral from the upstream region needs to be examined 
properly. 
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Upstream Contribution 


For purposes of calculation, the integration of equations 
(17) and (20) can be split into three regions: 

a) from upstream infinity to the beginning 
of the testsection (x=x s ) , 

b) from x s to the end of the testsection (x=x e ) , 

c) from x e to downstream infinity. 

Of these, as discussed above the contribution from region 
(c) is small for both test Mach number and angle of attack 
corrections, for practical size of the testsection lengths 
generally used. The region (b) is over which the wall 
pressure measurements are available, and is amenable for 
accurate calculation. However, the contribution from region 
(a) , for the angle of attack correction cannot be ignored and 
needs to be accounted properly. Since the wall pressure 
measurements are generally not available in region (a) , 
judicious interpolation across the front end or extrapolation 
of the range of measured pressure data to upstream infinity 
may be necessary in most of the interference calculation 
methods (ref. 24) . 

For large negative x, W 2 (x) tends to unity and the 
integral in equation (20) can be written as 



s? 

0 

O 

II 

1 

H 

H 

1 

H 

to 

(22) 

where 

*s 



*1 - /i( C p,u-Cp,l)/2] dx/h 

(23) 


— OC 

x s 

J 2 =/t( u m,u- u m,l>] dx/h 

(24) 


-oc 


The value of the integral I x depends on the method used 
to extrapolate the experimental data beyond x=x g . If the 
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difference between the measured upper and lower wall pressures 
at the most upstream location is small, the contribution from 
I-L can be ignored. However, when this is not possible, a 
suitable assumption for the variation of pressure coefficient 
with x has to be made. 

One of the methods suggested in reference (25) , is to use 
an exponential decay of the type Cp exp(x) . This was found 
to give satisfactory results for the test case experimental 
data obtained in the NAE perforated wall testsection on a BGK 
airfoil (ref. 1) . However, it may be noted that for this 
test case, the difference between the upper and lower wall 
pressure coefficients at the most extreme upstream location 
is 0.0001, and hence the calculated value of 1-^ is small in 
comparison with I 2 , and can be ignored. 

In the present investigation, an extrapolation based on 
the pressure distribution due to a vortex placed between two 
solid walls distance h apart has been considered. For this 
problem, the pressure coefficients on the top and bottom walls 
are given by 

C p (x,h/2)=-C p (x,-h/2) = - (7/0 h)/[ cosh (7rx//3h) ] (25) 

Using equation (25) , the difference between the top and 
bottom wall pressure coeficients can be written as 

C p ( x , h/ 2 ) - C p (x,-h/ 2 )=-( 27 / 0 h)/[cosh( 7 Tx// 3 h) ] (26) 

The equation (26) can be expressed in terms of the 
measured pressure coefficients at the most upstream location 
as 


C p (x , h/2) - C p (x ,-h/2) 


^p ( x s » h/ 2 ) ~ Cp(Xg, h/2) 


[cosh(7Tx s /j3h) ] 

(27) 

[ cosh (7rx/ 1 3 h) ] 
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For large negative x, equation (27) behaves as 


C p (x , h/2 ) - C p (x ,-h/2) 

= exp[ (x-x s ) (n/ph ) ] (28) 

C p (x s ,h/2) - C p (x s ,-h/2) 

It is expected that the extrapolation suggested by equa- 
tion (28) may represent an upper bound, since for ventilated 
walls, the actual value may be expected to be between the 
open jet and solid wall limits. The effect of ventilated 
walls is taken into account in an indirect manner by using 
the measured values of the pressure coefficients at x=x s in 
equation (28) . Substituting equation (28) in equation (23) , 
the integral 1-^ can be evaluated to give 

II = (P/210 <C piU -C p;1 ) x=xs (29) 

It may be observed that this value can be significant, if 
the wall pressure measurements are not extended far enough 
upstream so that the difference are small enough and can be 
ignored. Because of this, care should be exercised in deter- 
mining the location of x s for calculations, to ensure that the 
pressure measurements are not affected by any local flow 
conditions such as the beginning of ventilations. 

The integral I 2 represents contribution due to a vortex 
in free air, and forms a significant portion of the correc- 
tion. A closed form expression can be obtained for its 
value. For a vortex of strength 7 , the streamwise component 
of the induced velocity is given by 

u m (x,y) = [ (7/ 270 (py) ]/[x 2 +p 2 y 2 ] (30) 

and hence 

u m (x,h/2) - u m (x,-h/2) = (7/flj (|3h/2)/[x 2 +(j3h/2) 2 ] (31) 
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Substituting equation (31) in equation (24), the integral I 2 
can be evaluated to give 

I 2 = (cCi/2h7T) [ (7T/2) + tan -1 (2x s /(3h) ] (32) 

The variation of I 2 with x s /h is shown in figure (4) , and it 
may be seen that the contribution due to vortex singularity, 
from far upstream can form a significant portion of the cor- 
rection. For (c/h) = 0.25 and (x s /h)=l, the value of I 2 is 
about 0.8 degrees at a lift coefficient of about one. 


FOURWALL INTERFERENCE CORRECTION 

The procedure used for correcting the test data is the 
sequential approach suggested in reference (10) . First the 
test Mach number, the measured wall pressures and the model 
force coefficients are corrected for the sidewall boundary- 
layer effects. It may be noted that the sidewall boundary- 
layer correction theories presently used account only for the 
negative blockage caused by the enlargement of the streamtube 
due to favourable pressure gradient in the airfoil region. 
These sidewall boundary-layer corrected values are then used 
to calculate the top and bottom wall interference effects 
which results in an additional blockage correction, and a 
correction for the incidence due to lift interference. The 
measured force coefficients are then corrected for the change 
in dynamic head, using equation (19b) . 


DESCRIPTION OF THE COMPUTER PROGRAM 

The interference calculation method described above, has been 
incorporated into a Fortran computer program to facilitate 
calculation of wall interference effects on the airfoils 
tested in the Langley 0.3-m Transonic Cryogenic Tunnel. The 
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input details for the program are given in Appendix A. The 
first record of the input data determines one of the following 
options provided in the program: 

1) The required test data for calculating the wall inter- 
ference are read from the test data tapes. The array 
numbers where the values of the required test parameters 
are provided in a separate file. A typical example is 
shown in Appendix B. This file is incorporated into 
the main program by using the XEDIT facility. 

2) The required test parameters may be provided through the 
input namelist IDAT. 

3) The program can be run for the check case test data 
obtained on a BGK airfoil in the NAE perforated wall 
testsection. The required input data for this case 
are given in reference (1), and are incorporated in the 
present program. 

For the first two options, the various control parameters 
determining the range of integration, and the type of correc- 
tion to be applied are specified through the input namelist 
VALUS. The default values in the namelist have been set to 
correspond to cases generally used in evaluating interference 
for the airfoil test data in the 0.3-m TCT. 

It has been the experience with tests in the TCT that 
some of the wall pressures can be affected by leakage and 
hence need to be dropped to avoid erroneous calculation of the 
wall interference effects. For this purpose, provision has 
been made to skip up to five pressure points on both the top 
and bottom walls, through the parameters ISKIPT, ISKIPL, IT(5) 
and IL(5) . The wall pressures are interpolated linearly at 
the intermediate locations specified through the parameter 
XINCR. 
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The integrals in equations (10) and (11) required for the 
calculation of top and bottom wall interference effects are 
evaluated using trapezoidal rule. With the default options 
provided, the interference calculations are made for sidewall 
effects, top and bottom wall effects and combined fourwall ef- 
fects. However, if necessary, calculations can be made 
either for sidewall or top and bottom wall effects only by 
specifying the appropriate values for the parameters ISWL and 
ITB. 


The contribution to the angle of attack correction from 
upstream infinity to the beginning of the test-section due to 
extrapolation of the pressure coefficients is calculated using 
equation (29) . This can be suppressed by setting IEXTR=0. 

The calculation of the sidewall boundary-layer effects 
is straight forward and can be done by specifying test Mach 
number and the sidewall boundary-layer parameters DS (=25*/b) 
and the shape factor H. If a negative value is specified, 
the program calculates the required parameters using the 
empirical equations (8) and (9) . 

RESULTS AND DISCUSSION 
Calculations for the NAE Test Case 

For validation of the present program, first the top and 
bottom wall interference corrections were calculated for the 
test data given in reference (1) . This test data was obtained 
with a 10" chord BGK-1 airfoil model in the NAE 15"x 60" two- 
dimensional test section with perforated top and bottom walls. 
The test Mach number and Reynolds number were 0.784 and 21xl0 6 
respectively. The model was set at an incidence of 2.56 
degrees and the lift coefficient was 0.764. The results of 
interference corrections by various methods for this test case 
are summarised in reference (1) . 
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The required wall pressure data and other parameters for 
the NAE test case have been incorporated into the program and 
can be run by using option (3) discussed in Appendix (A) . 

The correction to the test Mach number obtained by the present 
method was -.015, and the correction to the angle of attack 
was -0.65 degrees. These values are compared in Table II, 
with the results of other methods taken from reference (1) . 

It may be noted that there is close agreement between the 
various methods and the present calculations. The correspond- 
ing corrections by the Kemp's method is -.017 for the Mach 
number and -.64 degrees for the angle of attack without 
upstream flow adjustment. With upstream flow adjustment, the 
correction to the angle of attack by the Kemp's method happens 
to be -0.89 degrees. 

As discussed earlier, the correction to the angle of 
attack is dependent on the type of extrapolation employed for 
extending the range of measured wall pressure data to upstream 
infinity, and the values of the top and bottom wall pressure 
coefficients at the most upstream location. For this test 
case, the difference between the top and bottom wall pressure 
coefficients at the most upstream location happens to be 
0.0001 and hence the contribution from the extrapolation as 
calculated from equation (29) happens to be about -.0006 
degrees and can be ignored. Hence the calculated angle of 
attack correction will not be much different, whether the 
wall pressures are extrapolated are not. In reference (25) , 
an exponential type of extrapolation for the wall pressures 
was used and the results were found to agree with the other 
methods. Perhaps, the small difference between the pressure 
coefficients at the most upstream location accounts for this 
agreement. 

However, the major upstream contribution to the angle of 
attack appears to come from the vortex singularity as calcu- 
lated from equation (32) . For a given testsection and 
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airfoil chord, this is directly proportional to the lift coef- 
ficent of the airfoil, and can be calculated independently 
of the wall pressure distribution. This happens to be a 
significant portion of the total correction. For the NAE 
test case, this value is -.26 degrees which forms about 40% 
of the total correction. This contribution appears to be a 
consequence of the finite length of the testsection. 

The foregoing discussion suggests that if the pressure 
measurements are available over a sufficiently long upstream 
distance so that the difference between the top and bottom 
wall pressure coefficents can be ignored in relation to other 
experimental uncertainties, all the methods give nearly same 
value of the correction for the both the Mach number and the 
angle of attack. When the difference between the most 
upstream top and wall pressure coefficients is significant, 
it may be necessary to consider extrapolation of the data. 
However, caution has to be exercised to ensure that the 
difference is really due to wall interference effect. Any 
extraneous local effect, or other experimental uncertainties 
can result in erroneous corrections for the angle of attack. 

Calculations for a Sample Case from 0.3-m TCT 

The sample case considered corresponds to the test data 
obtained on a 6” chord NACA 0012 airfoil in the Langley 0.3-m 
TCT. Results of fourwall interference calculations for this 
case using TWINTN4 have been given by Kemp in reference (10) . 
The test Mach number and Reynolds number were 0.701 and 3xl0 6 , 
respectively. The input data taken from reference (10) and 
used in the present calculations is given in Appendix C. The 
corrected Mach number obtained by the present combined four- 
wall calculations is 0.691 which is close to 0.689 obtained by 
Kemp's method. The corresponding values for the correction 
to angle of attack are -.25 and -.15 degrees for the present 
and Kemp's calculations respectively. It may be noted that 
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for this sample case, the pressure coefficients at the first 
measuring station located about four chords upstream are 
-0.005486 and -0.022372 respectively, on the top and bottom 
walls. Hence, the contribution from extrapolation of the 
pressure coefficients will be considerable and probably 
accounts for the difference in corrections obtained by the 
two methods. 

Calculations for Specific Test Programs in 0.3-m TCT 

The main objective of the present task was to automate 
the present fourwall interference correction procedure to 
facilitate a quick evaluation of the interference effects 
for all the data points in a specific test program. An 
assessment of these interference corrections can then form a 
basis for undertaking detailed interference calculations, if 
necessary, for specific data points likely to be affected by 
large transonic effects. 

The present program was validated by applying it to some 
of the airfoil test programs conducted in the 0.3-m TCT. The 
example chosen to demonstrate the procedure is the unpublished 
data by Mineck and Lawing, on a 6" chord, 12% thick symmetric 
super-critical airfoil model in the 8"x 24" slotted wall 
testsection. The sequence of control cards and the input 
data, and typical output result obtained are shown in Appendix 
D. This particular example of a symmetric airfoil was chosen 
to identify the problems associated in making interference 
corrections and the associated uncertainties. As described 
in Appendix A, the program is run in option (1) , by specify- 
ing the test and run numbers only. The program calculates 
the interference corrections for all the test points in the 
specified run number, and the results are printed in the out- 
put format shown in Appendix D. The corrections to the Mach 
number and the angle of attack (in degrees), and the corrected 
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values of the test parameters are printed for each type of 
correction: i.e., sidewall, top and bottom wall, and combined 
fourwall interference. 

Correction for Blockage 

The correction to the test section Mach number for a 
typical run in the 0.3-m TCT is shown in figure (5), at a Mach 
number of 0.6 and for various lift coefficients. It may be 
noted that the correction for the sidewall boundary-layer 
effects as obtained by equation (4), for the 0.3-m TCT is 
significant amounting to about -0.014. The correction due to 
top and bottom wall effects is about 0.002 and does not seem 
to vary significantly over the range of lift coefficients from 
-0.5 to 0.5. The relatively small blockage correction due to 
top and bottom wall effects is expected for the 0.3-m TCT 
slotted wall testsection which has been designed for low 
blockage effects using the method of reference (21) . It may 
be noted that such slotted wall tunnels designed for low 
blockage effects, may introduce significant correction to the 
angle of attack. However, considering the uncertainties 
involved in determining the true angle of attack in two- 
dimensional airfoil testing, designing for low blockage 
effects is an attractive feature. 

Correction for Lift Interference 

It has to be noted that the calculated values of the 
correction for the angle of attack needs to be accounted for 
any testsection flow inclination. It is desirable that this 
information is to be obtained from empty tunnel calibration. 
However, for a symmetrical airfoil this can be deduced from 
lift coefficient versus angle of attack curves. For the 
symmetrical airfoil tested in the 0.3-m TCT, these are shown 
plotted in figure (6) for Mach numbers 0.5, 0.6, 0.76 and 0.8. 
From these lift curve data, it appears that an empty test- 
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section flow inclination of about -0.1 degree is required for 
the test data. 

Again, for a symmetrical airfoil at zero lift, the top 
and bottom wall pressure signatures has to be identical and 
hence the calculated correction for the angle of attack has to 
be zero. However, any local flow conditions at the walls or 
other effects may yield a non-zero correction, which can be 
considered as a tare correction due to measured wall pressures 
at zero lift. This is demonstrated in figure (7) by plotting 
the calculated correction for the angle of attack for various 
lift coefficients at a Mach number of 0.6. The difference 
between the corrections obtained by extrapolating the measured 
pressures to upstream infinity, and without extrapolation is 
not significant for the case considered. The calculated 
value of the correction varies nearly linearly over the range 
of lift coefficients from -0.5 to 0.5. The tare correction 
due to measured wall pressures at zero lift happens to be 
about 0.1 degree. Assuming this tare correction remains 
constant with lift, this value needs to be substracted from 
the calculated value of correction using the measured pressure 
distribution. An example of the application of this correc- 
tion procedure is illustrated in figure (8) , corresponding to 
a test Mach number of 0.6. The corrected lift curve is 
closer to empirically correlated Davis-Moore theory (ref. 21 
and ref. 26) for positive lift coefficients. For negative 
lift coefficients the difference between the two methods is 
noticeable. 


Application to a Cambered Airfoil 

The present correction procedure was applied to a recent 
test on a cambered super-critical airfoil in the 0.3-m TCT 
slotted wall test section. Detailed wall interference cal- 
culations by using the TWINTN4 code have been made recently+ 
for this airfoil, and hence it was thought that a comparison 
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of the present method with the calculations of TWINTN4 would 
provide a better assessment of the method. in figures (9) 
and (10), typical results obtained by the present method are 
shown for two Mach numbers, 0.6 and 0.73 respectively, for a 
chord Reynolds number of 30 million. Good agreement with the 
TWINTN4 results was found for these cases and for many other 
test conditions. This suggests that the present method can 
be employed to get first order interference corrections for a 
typical test program, which can subsequently be used to deter- 
mine test conditions requiring more detailed evaluation using 
TWINTN4. With the present computer program, the calculations 
can be done for all the test points with little computational 
effort. 


Practical Problems 

From the several examples considered above, and the test case 
data shown in Table II, it follows that for a given subcriti- 
cal wall pressure distribution, the interference corrections 
obtained by the present program and the various other methods 
are nearly same. The application of most of these methods is 
relatively straightforward except for the method of Kemp which 
solves the non-linear transonic problem. However, when 
applied to a specific test program, several difficulties can 
arise, mainly due to experimental limitations. This leads to 
uncertainties in the calculated value of the interference cor- 
rections and often it is difficult to overcome these problems 
by making refined calculations. Some of these problems have 
been addressed by Smith (ref. 24) . With particular reference 
to measurements in the 0.3-m TCT, these problems are mainly 

a) The upstream wall pressure measurements, depending on the 
location of the presssure orifices, can be affected by the 
local flow conditions in the vicinity of the beginning of the 
top and bottom wall ventilations. This can introduce spurious 
pressure signatures not related to the model perturbations. 
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Often judgement has to be exercised on deciding the most 
probable location where the measurements are not affected by 
the local flow conditions. 

b) Ideally, it is desirable that the pressure measurements 
are made at a distance away from the wall by using pressure 
rails. However, for practical reasons, the measurements are 
often made on the slats. The accuracy of using the slat 
measurements instead of the pressure rail measurements has 
been examined by Smith (ref. 24) by making both the measure- 
ments in the NLR Pilot Tunnel. For the CAST-7 airfoil with a 
chord to tunnel height ratio of about 0.33, it was found that 
the differences between the two measurements were not sig- 
nificant, and the calculated values of the corrections were 
well within the overall experimental accuracies. However, 
this is a factor which depends largely on the ratio of slot 
spacing to tunnel height, and it is difficult to draw any firm 
conclusion in the absence of experimental data for a specific 
facility. 

c) Due to limited testsection length, extrapolation or 

interpolation of the measured wall pressure data is necessary 
for interference calculation methods requiring boundary data. 
This is particularly important for determining the correction 
for the angle of attack. The suggested method of using the 

wall pressure distribution due to a vortex between solid walls 
provides a theoretical basis for extrapolation and can be 
expected to give an upper bound, since for an open jet the 
pressures will be identically same all along the boundary. 
However, it is desirable to check the angle of attack correc- 
tions obtained with and without extrapolation of the pressures 
and ensure that there is not significant difference between 
the two values. The present method provides a means for 
making such quick checks and sensitivity analysis of the 
interference corrections. 
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d) Far upstream, the difference between the wall pressures 
and the freestream static pressure tends to become small. 

Since the correction for the angle of attack depends on the 
difference between the top and bottom wall pressures, care 
needs to be taken to measure these differences accurately. 

If there is considerable noise, the measured data may have to 
be smoothened appropriately. 

e) Ideally, for a symmetrical airfoil, one would expect the 
correction for the angle of attack to be of the same magnitude 
but of opposite sign when the airfoil is at negative angle of 
attack. However, it may be difficult to achieve this in 
ventilated walls due changing local flow conditions. 

f) In ventilated wall tunnels, changes in upstream flow 
inclination can occur with change in lift coefficient. This 
problem has been addressed in detail by Kemp. It appears 
that the scope of the present program can be extended, if 

a suitable correction based on the upstream pressure 
measurements can be made to account for the upstream flow 
inclination. 


CONCLUDING REMARKS 

A simplified fourwall interference calculation procedure 
has been developed to correct the airfoil data obtained in the 
0.3-m TCT. The procedure was applied to typical test cases 
and airfoil test programs in the TCT, and good agreeement was 
observed with the results of other methods. While the 
application of the various interference calculation methods 
is straightforward, it appears that the practical limitations 
and uncertainties associated with the experimental data can 
impose limits on the accuracy to which the interference 
corrections can be assessed. 
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The correction to the test Mach number due to blockage 
effects is not much affected by the uncertainties associated 
with extending the range of the wall pressure distribution 
beyond the measured limits. However, the same argument does 
not hold true for the angle of attack correction in ventilated 
wall tunnels. This often forces the tunnel engineer to dis- 
regard the experimental angle of attack. To quote from 
reference (27) , 

-"For two-dimensional tests the wall ventilation is 
generally configured to minimize blockage interference which 
predominantly affects the Mach, number of the freestream. 

The lift interference effects are so large that the experi- 
mental angle of attack is disregarded and section normal force 
coefficient is emphasized". 

For the 0.3-m TCT, significant portion of the blockage 
correction comes from the sidewall boundary-layer effects, 
since the top and bottom slotted walls are designed for low 
blockage effects. Experience with the analysis of the 0.3-m 
TCT airfoil data (refs. 28, 29) suggests that a correction for 
the test Mach number based on the one-dimensional effect of 
sidewall boundary-layer effects is often adequate to give 
satisfactory results, for the 6" chord models generally tested 
in the TCT. 

It appears that the empirically correlated Davis-Moore 
theory often gives acceptable correction for the angle of 
attack (e.g., fig. 8) for the 0.3-m TCT airfoil data. 

However, as has been observed in reference (26) , this may not 
be true in general, and methods based on measured boundary 
conditions are superior and more reliable. However, while 
arriving at corrections for the angle of attack, care needs to 
be taken in ascertaining the quality of the wall pressure 
data. 
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The present method uses simple singularities for model 
representation, and the trapezoidal rule for integration of 
wall pressures. This was done to keep the calculation 
simple, so that, the method can be adopted for making on-line 
corrections for the airfoil tests in the TCT. However, the 
method can be improved by using a better model representation, 
and integrating the wall pressure using a spline fit for the 
measured data. Further, the scope of the approach can be 
improved if upstream flow inclination can be properly ac- 
counted . 

Good agreement of the present method with the results of 
TWINTN4 code suggests that the method can be used to obtain a 
quick correction for all the data points in an airfoil test 
program. A detailed evaluation of the interference may then 
be made for selected conditions using TWINTN4 . 

The present approach can be used directly for calculating 
the correction for the angle of attack, even with upstream 
sidewall boundary-layer removal since it introduces identical 
pressure signatures on both the top and bottom walls. How- 
ever, for the Mach number, a tare correction to account for 
mass removal effects will be required. This can be obtained 
either from empty tunnel calibration with sidewall boundary- 
layer removal or from theoretical considerations. 
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Table I 


Summary of 2-D Wall Correction Methods 


Method 


Input 

Boundary 

Model 

Analysis 

Mokry 

(Ref .13) 

Pressures 

Forces 

Subsonic 

Capalier et al. 
(Ref. 12) 

Pressures 

Forces 

Subsonic 

Ashill 

(Refs. 

& Weeks 
16 & 17) 

Pressures & 
Flow Incln. 

None 

Subsonic 

Sawada 

(Refs. 14 & 15) 

Pressures 

Forces 

Subsonic 

Smith 

(Panel 

(Ref. 4) 
method) 

Pressures 

Pressures, 
Wake drag 

Subsonic 

Kemp 

(Finite 

(Ref. 3) 
difference) 

Pressures 

Pressures 

Transonic 

Murman 

(Finite 

(Ref. 2) 
Difference) 

Pressures 

Pressures 

Transonic 
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Table II 


Comparison of Corrections* for the NAE Test Case 


Correction for 


Method 

Mach No. 

Incidence 

Mokry, Ohman (Ref. 13) 

-.015 

-.67 deg. 

Capalier et al. (Ref. 12) 

-.015 

-.67 deg 

Gopinath (Ref. 25) 

-.017 

-.67 deg. 

Smith (Ref. 4) 

-.015 

-.59 deg. 

Sawada (Ref. 14, 15) 

— 

-.58 deg. 

Kemp (Ref. 11) 

-.017 

-.64 deg. 



-.89 deg. 

Present Calculations 

-.015 

-.65 deg. 


* taken from reference (1) 

+ with upstream flow angle adjustment 
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Figure 4 : Contribution 


from 



x s/ h 

to angle of attack due to vortex 
: infinity to x s (See equation 32) . 
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Figure 7 : Calculated correction for the angle of attack for 

the symmetrical supercritical airfoil. 
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Figure 8: Corrected and measured lift coefficients for the 

symmetrical supercritical airfoil. 
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Application of the present method to a cambered 
airfoil test (subcritical case) . 
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APPENDIX A: PROGRAM INPUT DETAILS 


Option 1: Calculations for a given Test. 


Namelist VALUS 

Record 1. 00 (Columns 1-2); Input data for wall pressures 

and airfoil force/pressure 
coefficients read from test 
data files. 

Record 2 Test number ( 13) ( Run N umbe r ( I 2 ) , D S ( F 1 0 . 6 ) , S H ( F 1 0 . 6 ) 

Record 3 *-,,** 

etc., 

Last 00000 (Columns 1*5), DS, SH (2F10.6) 

/ E 0 R 


Option 2: Calculations for a given case. 


Namelist 

VALUS 



Record 

1 . 

01 

(Col umns 1 

*2) 

Record 

2 . 

TITLE 

(Col umns 1 * 

80) 

Record 
/ E 0 R 

3 . 

$ I D A T 

(From Col umn 

2) 


Record 
/ E 0 R 


Option 3: Calculations for Check Case (NAE Data) 


02 (Col umns 1*2) 
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Parameters in Namelist VALUS 


XSTART : 

XEND : 

X I NCR : 

AREA : 
WIDTH : 
H : 

CHORD : 
XL : 

I AR : 

WVL : 

I SWL : 


IEXTR : 


I SK I PT 
I SK I PB 
I T ( 5 ) 

I B ( 5 ) 

I REND 
I PEND 
I TB 


(default values of the parameters are shown in paranthesis) 

Upstream value of x for integration of wall 
pressures (-24.5) 

Downstream value of x for integration of wall 
pressures ( 23.5) 

Value of increment for X interpoalting the inter- 
mediate values ( 2.0) 

Cross-Sectional area of airfoil ( 3.0 sq" ) 

Width of the tunnel ( 8.0") 

Distance between the top and bottom walls ( 24 " ) 

Airfoil Chord (6.0") 

Distance of airfoil leading edge from turntable 
center ( 2.16") 

1, for including aspect ratio effects (Default) 

0 , aspect ratio effects not included 

Length scale representing airfoil chord. Required 

if aspect ratio correction is made. (2) 

=0, No sidewall boundary-layer correction. 

=1, Murthy's Correction method. (Default option) 

= 2, B a r nw e l l - S e wa l l correction method. 

=1, Upstream most top and bottom wall pressures are 
extrapolated to infinity. (Default option) 

=0, No extrapolation of wall pressures. 

No. of upper wall ports to be skipped (Max. 5) 

No. of lower wall ports to be skipped (Max. 5) 

Array, Port numbers to be skipped on top wall (0) 

Array, Port numbers to be skipped on bottom wall (0) 

Last Run number on the data files (100) 

Last Point number on thedata files (100) 

=0, No top/bottom wall corrections applied 

*0, top/bottom wall correction applied (default) 
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Parameters in Namelist IDAT 


EM 

Test 

Mach number 




ALPHD 

Angle of Attack 

(Deg . ) 



CL 

Lift 

Coefficient 




CD 

Drag 

Coefficient 




DS 

Sidewall boundary - layer thickness (2<$*/b) 


SH 

Sidewall boundary - layer shape factor 


RE 

Chord Reynolds number 



NUA F 

No. 

of orifices 

on airfoil 

upper s u r f ac e ( M ax . 60 ) 

NLAF 

No. 

of orifices 

on airfoil 

lower surf ace(Max.60) 

NU 

No . 

of orifices 

on the top 

wall (Max. 

60) 

NL 

No . 

of orifices 

on the bottom wall (Max. 

60) 

XAFUC60) 

x/c 

location of 

orifices on 

airfoil upper 

surface 

XAFLC60) 

x/c 

location of 

orifices on 

airfoil lower 

surface 

X WU ( 6 0 ) 

x/c 

location of 

orifices on 

top wall 


X W L ( 6 0 ) 

x/c 

location of 

orifices on 

bottom wall 


CPAFUC60) 

Cp 

values on ai 

rfoil upper 

surface 


CPAFLC60) 

Cp 

values on airfoil lower 

surface 


CPWUC60) 

Cp 

values on top wall 



CPWL ( 60 ) 

Cp 

values on bottom wall 




* 

Not used in the program. 
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Appendix B 

Example of file for Array Numbers for Test Parameters 


ARRAY NUMBERS RELEVANT TO THE TEST ARE INSERTED HERE 
ARRAY NUMBERS FOR TEST190 (ASH) 

N P ( 1)= 1076 

N P ( 2 ) = 1078 
N P ( 3 ) = 1079 
N P ( 4) = 1065 
N P ( 5 ) = 1220 
N P ( 6 ) = 1066 
N P ( 7 ) = 1067 

NP ( 1 2 )= 45 

N0U.N0L: NO. OF ORIFICES ON AIRFOIL UPPER, LOWER SURFACES 

NOT, NOB: NO. OF ORIFICES ON TUNNEL TOP AND BOTTOM WALLS 


c 

I STX 

, I CPST : 

STARTING 

ARRAY 

NUMBER 

FOR 

X 

AND 

CP, 

(TOP 

WALL ) 

c 

I SBX 

, I CPSB : 

START I NG 

ARRAY 

NUMBER 

FOR 

X 

AND 

CP, 

(TOP 

WALL ) 

c 

I USX 

, I CPUS : 

START I NG 

ARRAY 

NUMBER 

FOR 

X 

AND 

CP, 

(AIRFOIL 

US ) 

c 

p 

I LSX 

, ICPLS: 

STARTINF 

ARRAY 

NUMBER 

FOR 

X 

AND 

CP, 

(AIRFOIL 

LS) 

L 

N 0 U = 

25 S NOL 

=28 $ NOT 

= 26 $ 

N 0 B = 2 6 









ISTX = 1829 $ I C P S T = 570 
ISBX = 1566 $ I C P S B = 596 
IUSX=1983 $ ICPUS= 492 
I LSX = 1350 $ I C P L S = 517 
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Appendix C 

Input Data for the Sample Case (Ref. NASA CR-3777) 


$ V A L U S XSTART=-26.5, XEND = 23.5, X I N C R = 2 . 0 , I A R = 0 , I S WL = 2 , $ E N D 

01 

DATA FROM NASA CR-3777, P-24, TUINTN4: REFrWILLIAM B KEMP, JR. 

$IDAT EM=. 70 1,CL=. 2204, CD = . 0076, DS = 0. 01543, SH = 1. 5042, RE=6,ALPHD = 0.0, 

N U = 2 6 , N L = 2 8 , NUAF = 24, N L A F = 2 4 , 

CPUU( 1 )= - .54861E-02, 


• . 1 2283E - 02 

f - . 98 1 2 2 E - 0 3 , 


.62431E-03, 

- .52061E-03, 

.86244E-03, 

.11078E-03 

. 13516E-02 

, - .37153E-02, 

- 

. 84227E - 02, 

- . 19627E - 01 , 

- .36989E-01 , 

- . 4 8 0 8 2 E - 0 1 

- .60229E-01 

, - . 5 92 75 E - 0 1 , 

- 

. 5 0005 E - 01 , 

- .50088E - 01 , 

- . 321 47E-01 , 

- .33204E-01 

• . 2 8 0 5 6 E - 0 1 

, - .32320E-01 , 

- 

. 3 6 9 63 E - 0 1 , 

- .371 1 7E - 01 , 

- . 36747E - 0 1 , 

- . 3 9867E - 0 1 

- . 3 9 5 5 2 E - 0 1 







C P W L ( 1 )= - 

.22372E-01 , 






. 69 1 73 E - 0 3 

, * .35306E-02, 


.47236E-02, 

. 14234E - 01 , 

. 98306E - 02 , 

. 13424E-01 

. 22097E - 01 

, .18123E-01, 


. 22969E - 01 , 

. 9 78 8 8 E - 02, 

- . 56598E - 02, 

.21510E-02 

. 2 0 2 83 E - 0 3 

, .14175E-02, 


.84934E-03, 

- . 8 8 3 0 7 E - 03, 

- .23261E-02, 

- .70418E-02 

- .92554E-02 

, -.15335E-01, 

- 

. 1 6095E - 01 , 

-.18531E-01, 

- . 1 8062E - 01 , 

- . 18691E-01 

• . 18373E-01 

, -.17930E-01, 

- 

. 4 5 1 4 2 E - 0 1 , 




XWU( 1 )= -26.5, -24.5, -22 

. 

5, -20.5, - 

18.5, -16.5, 

-14.5, -12.5, 


- 

10.5, -8.5, 

6 

.5, *4.5, 

-2.5, -0.5, 

1.5, 3.5, 



5.5, 7.5, 

21.5, 23.5, 

9 

.5, 11.5, 

13.5, 15.5, 

17.5, 19.5, 


X W L ( 1 )= - 26.5 , - 24.5 , - 22 

. 

5, -20.5, - 

18.5, -16.5, 

-14.5, -12.5, 



10.5, -8.5, - 

6 

.5, -4.5, 

-2.5, -0.5, 

1.5, 3.5, 



5.5, 7.5, 9 

. 

88, 11.88, 

13.88, 15.88, 

17.88, 19.88, 


CM 

CO 

oo 

CM 

CM 

OO 

CO 

O 

CM 

. 

88, 26.88, 

28.88, 



XAFU( 1 ) = 

.22477E- 03, 






. 14361E-01 

, . 2 8 2 2 1 E - 0 1 , 


. 5 2 9 70 E - 0 1 , 

.77977E-01 , 

. 1 0296E + 00, 

. 15231E+0Q 

.20155E+00 

, -25228E+00, 


.30200E+00, 

■35197E+00, 

.40145E+00, 

.45146E+00 

.50149E+00 

, .55082E+00, 


.60089E+00, 

.65075E+00, 

.70067E+00, 

.74982E+00 

. 80091 E + 00 

, -85039E+00, 


.89945E+00, 

.94878E+00, 

.987503 , 


CPAFUC 1 ) = 

. 10741E + 01 , 






• . 32689E + 00 

, - -61743E + 00, 

- 

.83530E+00, 

- .92300E + 00, 

- . 97980E + 00 , 

- .94881 E+00 

- . 84 726E + 00 

, -.70543E+00, 

• 

.60428E+00, 

- .53857E + 00, 

- .48697E + 00, 

- .42540E+00 

- .37318E + 00 

, - .31092E+00, 

- 

.26626E+00, 

- . 22483E + 00, 

- . 1 7602E + 00, 

- . 1 2971 E+00 

- . 74 0 6 2 E - 0 1 

, • .20218E-01, 


.4471 8E-01 , 

.12041E+00, 

.21364 , 


XAFL< 1 ) = 

.22477E-03, 






. 1 1 7 1 8 E - 0 1 

, . 243 4 1 E - 0 1 , 


. 4 993 6 E - 01 , 

.74166E-01 , 

.98549E-01 , 

. 14860E+00 

. 1 9846E + 00 

, -24811E+00, 


.29807E+00, 

.34861 E+00, 

•39795E+00, 

.44826E+00 

.49806E+00 

, . 54869E+00, 


.59798E+00, 

.64853E+00, 

.69885E+00, 

. 74836E+00 

.79859E+00 

, .84924 E+00, 


. 899 1 4 E + 00 , 

■94766E+Q0, 

.987503 , 


C P A F L ( 1 ) = 

. 1 0741 E + 01 , 






.49190E+00 

, . 17864E + 00, 

• 

.10019E+00, 

- .21237E + 00, 

- .29413E + 00, 

- .36746E+00 

- .38417E+00 

, - .3 73 48E + 00, 

• 

.34969E+00, 

-.32701E+00, 

- . 301 94E + 00, 

- . 26462E+00 

- .23905E+00 

, --20599E+00, 

- 

.17514E+00, 

- .15402E + 00, 

- . 1 1 758E + 00 , 

- .85763E-01 

- . 4 0 8 1 4 E - 0 1 

, - .83148E-02, 


. 2 9 5 8 5 E - 0 1 , 

.1 1052E+00, 

.255044 , 



SEND 
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« .r 


« * * n 


Results for the Sample Case 


DATA FROM NASA CR-3777, P-24, TWINTN4: 


REF:WILLIAM B KEMP, 


BARNUELL/SEWALL CORRN (DOES NOT ACCOUNT FOR AR) 


DELTA M 

• .0139 

D 

ALPHA 

0 .0000 

M 

CORR 

.6871 

CL 

CORR 

. 2234 

CD 

CORR 

.007702 

AL 

CORR 

0 .0000 

TOP - BOT 

WALL CORRE 

DELTA M 

. 004 1 

D 

ALPHA 

• . 2588 

M 

CORR 

. 705 1 

CL 

CORR 

.2186 

CD 

CORR 

.007539 

AL 

CORR 

• . 2588 


BOTH SIDEWALL & TOP-BOT WALL CORRECTION 


DELTA M 

.0042 

D 

ALPHA 

- .2516 

M 

CORR 

.6913 

CL 

CORR 

.2215 

CD 

CORR 

.007637 

AL 

CORR 

- .2516 


JR . 
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Appendix D 


Control Cards and Input for a typical 0.3-m TCT Test Program 


FOURADX # T 1 000 , CM170000. 

USER , USERNUM, PASSWOR . 

CHARGE , XXXXXX , XXX . 

GET f FETCH/UN=474750C. 

GET, J190001/UN=USERNUM. 

FETCH(TAPEI) 

GET # FOURADX f ARAY 1 90 . 

XE # FOURADX . 

F T N , I=FOURADX,L=0,PL=30Q0. 

ATTACH(FTNMLIB/UN=LIBRARY,NA) 

LDSET(LIB=FTNMLIB # PRESETA=INDEF) 

L GO . 

REPLACE, TAPE10=OUPTFLN. 

DELIVER .XXXXXXX XXXXXX 

GET,OUPTFLN- 

COPY, OUPTFLN, OUTPUT. 

EXIT. 

REPLACE, TAPE10=ERRADX. 

GET , ERRADX . 

COPY, ERRADX, OUTPUT. 

/ E 0 R 
J 1 9000 1 
/ E O R 
DEL; 

L/IDAT=00/;READ ARAYl90;T;D/%/*;E 
/ E O R 

SVALUS I T ( 1 > = 10, ISKIPT = 1 , 1SKIPL = 1 , I L( 1 ) = 1 , I T B = 1 , IAR = 0, IREND = 1 , IPEND = 12,$END 

00 

19001 -.0206 1.4448 

00000 .0000 .0000 

/EOF 
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Typical output for a 0.3-m TCT Test 


$ V A L U S 

XSTART = -.245E+02, 

XINCR = . 2 E + 0 1 , 

XEND = . 23 5 E + 0 2 , 

AREA = . 3 E + 0 1 , 

H = . 2 4 E + 0 2 , 

XL = .216E+01, 

CHORD = . 6 E + 0 1 , 

I SK I PT = 1 , 

I SK I PL = 1, 

WIDTH = . 8 E + 0 1 , 

I AR = 0 , 

WVL a . 2E + 01 , 

I SWL = 1 , 

I TB = 1 , 

IUAF = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

I L A F = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

ISKPAFU = 15, 

I SKPAFL = 15, 

DS * . 2E - 01 , 

SH = . 14E + 01 , 

I REND = 1 , 

IPEND = 12, 

IT =10,0,0,0,0, 

IL = 1 , 0, 0, 0, 0, 

I EXTR = 1 , 

SEND 


0 , 

0 . 


contd/ 
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TEST 

190. 

190. 

190. 

190. 

190. 

190. 

RUN 

1 . 

1 . 

1 . 

1 . 

1 . 

1 . 

POINT 

1 . 

2. 

3. 

4. 

5 . 

6. 




MEASURED 

VALUES 



ALPHAD 

.0102 

1.0122 

2 . 0060 

2.0060 

4.0120 

. 0000 

MACH NO. 

.6018 

.6013 

. 6006 

.6005 

.6026 

.6015 

REX 1 0 - 6 

5.9873 

5.9956 

5 . 9882 

5.9863 

5.9934 

5.9856 

CL 

- .0260 

. 0926 

.1857 

. 2005 

. 4425 

- . 0197 

CD 1 

.007225 

.007197 

. 007345 

.007331 

. 009376 

.007220 



SIDEWALL 

BL PARAMETERS (CALCULATED) 


2DS/B 

. 0207 

. 0207 

. 0207 

. 0207 

. 0207 

. 0207 

H 

1 .4457 

1 . 4457 

1 . 4457 

1 .4457 

1 .4457 

1.4457 



S I DEWALL 

BL CORRECTION (HURTHY) 


DELTA M 

- . 0140 

- .0140 

- . 0140 

- .0140 

- .0140 

- . 0140 

D ALPHA 

0.0000 

0.0000 

0 .0000 

0 .0000 

0 . 0000 

0 .0000 

M CORR 

.5878 

.5873 

. 5867 

. 5866 

. 5886 

. 5876 

CL CORR 

- . 0266 

. 0948 

.1901 

. 2053 

.4530 

- . 020 1 

CD CORR 

.007397 

. 007368 

. 00751 9 

.007505 

.009598 

.007391 

AL CORR 

.0102 

1.0122 

2.0060 

2.0060 

4.0120 

. 0000 



TOP - BOT 

WALL CORRECTION ONLY 



DELTA M 

. 003 1 

. 0043 

. 0022 

. 0027 

. 0032 

. 0022 

D ALPHA 

. 0957 

. 0955 

- . 0684 

- . 0825 

- . 3253 

. 1308 

M CORR 

. 6049 

. 6055 

. 6028 

. 6032 

. 6058 

. 6037 

CL CORR 

- . 0258 

.0917 

.1846 

. 1 991 

.4390 

- .0195 

CD CORR 

.007169 

. 0071 20 

.007304 

.007281 

. 00930 1 

.007180 

AL CORR 

. 1058 

1 . 1077 

1 .9376 

1 .9235 

3 . 6868 

. 1308 


BOTH 

SIDEWALL 

& TOP-BOT 

WALL CORRECTION 


DELTA M 

.0030 

. 0042 

. 0022 

. 0027 

. 003 1 

. 002 1 

D ALPHA 

. 0983 

. 0888 

- .0835 

- . 0989 

-.3618 

. 1333 

M CORR 

.5909 

.5915 

.5888 

. 5892 

.5917 

. 5897 

CL CORR 

- .0264 

.0938 

.1890 

.2038 

.4493 

- . 0200 

CD CORR 

.007339 

. 007288 

.007477 

. 007453 

. 009520 

.007350 

AL CORR 

.1085 

1 . 1009 

1 .9225 

1 .9071 

3.6503 

.1333 







c on t d . / 
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TEST 

190 . 

190. 

190. 

190. 

RUN 

1 . 

1 . 

1 . 

1 . 

POINT 

7 . 

10. 

1 1 . 

12. 




MEASURED 

VALUES 

ALPHAD 

• 1 .0081 

- 1 . 9958 

-3.9917 

. 0000 

MACH NO. 

.6011 

. 6003 

.6010 

. 5997 

REX10- 6 

5 .9826 

5.9917 

5.9997 

5 . 9923 

CL 

• . 1351 

• . 2489 

- .4808 

- . 0228 

CD 1 

.007425 

. 007872 

. 009376 

.007370 



SIDEWALL 

BL PARAMETERS (CALCULATED) 

2DS/B 

. 0207 

. 0207 

. 0207 

.0207 

H 

1 . 4457 

1 .4457 

1.4457 

1 .4457 



S I DEWALL 

BL CORRECTION (MURTHY) 

DELTA M 

- . 0140 

* .0140 

- .0140 

- .0139 

D ALPHA 

0 . 0000 

0.0000 

0 .0000 

0.0000 

M CORR 

. 5871 

. 5863 

. 5870 

. 5857 

CL CORR 

* .1384 

- . 2548 

- . 4922 

- .0234 

CD CORR 

.007601 

. 008059 

.009599 

.007545 

AL CORR 

• 1 . 0081 

- 1 . 9958 

-3.9917 

. 0000 



TOP - BO T 

WALL CORRECTION ONLY 

DELTA M 

. 0026 

. 0025 

. 0023 

. 0026 

D ALPHA 

.1628 

.2516 

. 4498 

.0917 

M CORR 

.6037 

. 6027 

. 6033 

. 6023 

CL CORR 

- .1342 

- .2473 

- .4780 

- . 0227 

CD CORR 

. 0073 75 

. 007823 

.00932 1 

.007321 

AL CORR 

- . 8453 

- 1 . 7442 

-3.5419 

.0917 


BOTH 

SIDEWALL 

& TOP-BOT 

WALL CORRECTION 

DELTA M 

. 0026 

. 0024 

. 0023 

. 0026 

D ALPHA 

. 1 744 

.2726 

.4904 

.0939 

M CORR 

. 5897 

. 5887 

. 5893 

. 5883 

CL CORR 

- .1374 

- . 2532 

- .4892 

- . 0232 

CD CORR 

.007550 

. 008008 

.009541 

.007494 

AL CORR 

• . 8337 

- 1 . 7233 

-3.5013 

. 0939 
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